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Abstract 

We calculate the ratio of shear viscosity to entropy density for charged black 
branes in third order Lovelock theory. For chargeless black branes, the result turns 
out to be consistent with the prediction made in l arXivi 0808.3498[hep — th]. We 
find that, the third order Lovelock gravity term does not contribute to causality 
violation unlike the Gauss-Bonnet term. The stability of the black brane again 
requires the value of the Lovelock coupling constant to be bounded by 1/4 in the 
infinite dimensionality limit. 

1 Introduction 

The AdS/CFT correspondence [1-3] provides an interesting theoretical framework for 
studying relativistic hydrodynamics of strongly coupled gauge theories. The result of 
RHIC experiment on the viscosity /entropy ratio turns out to be in favor of the predic- 
tion of AdS/CFT [4-6]. Some attempt has been made to map the entire process of RHIC 
experiment in terms of gravity dual [7]. The way to include chemical potential in the 
theory was figured out in [8,9]. 
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It had been conjectured that the viscosity value of theories with gravity dual may give 
a lower bound for the 'h/s = for all possible liquid [14]. However, in the presence of 
higher- derivative gravity corrections, the viscosity bound and causality are also violated 
as a consequence [15-18]. The ratio of shear viscosity to entropy density are of particular 
interest in higher derivative gravity duals because those higher derivative terms can be 
regarded as generated from stringy corrections given the vastness of the string landscape. 
In [19-21], the authors computed t]/ s for general gravity duals by determining the ratio of 
two effective gravitational couplings. The r]/ s in presence of arbitrary B? and B? terms in 
the bulk action were calculated in [22]. 

The higher derivative terms may be a source of inconsistencies because higher powers 
of curvature could give rise to fourth or even sixth order differential equation for the 
metric, and in general would introduce ghosts and violate unitarity. Zwiebach and Zumino 
[23, 24] found that ghosts can be avoided if the higher derivative terms only consist of the 
dimensional continuations of the Euler densities, leading to second order field equations 
for the metric. These theories are the so called Lovelock gravity [25]. The zeroth order 
of Lovelock gravity correspondences to the cosmological constant. The first order is the 
Einstein equation and the second order correspondences to Gauss-Bonnet theory. Higher 
derivative effects on rj/s in the presence of a chemical potential have been discussed in 
[26-33]. In this paper, we discuss shear viscosity in third order Lovelock gravity. 

Our motivation for this paper is based on the following facts: 

1) . Although people expect that rj/s might receive corrections from third and higher 
order Lovelock terms, it was conjectured in [19] that rj/s gets no corrections at all for higher 
order Lovelock terms except the Gauss-Bonnet terms. In this paper, we compute rj/s for 
third order Lovelock gravity directly by using the standard method developed in [4, 5] and 
compare our result with that of [19]. 

2) . In [15] and [16], the authors showed that if we consider the Gauss-Bonnet correction 
to Einstein equation, the viscosity bound is violated in the hydrodynamics regime. More- 
over, causality violation happens in the high frequency regime {k^ oo), which implies 
that theories in that regime are pathological [It]*. In [26] and [27], some of us considered 
medium effect and the higher derivative correction simultaneously by adding charge and 
Gauss-Bonnet terms and found that the viscosity bound as well as causality violation is 



*The causality issue in Gauss-Bonnet gravity was further studied in [34]. 
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not changed by the charge. After adding the third order Lovelock terms, the causahty 
structure would be different from that of the Gauss-Bonnet gravity. 

3). It is worth to studying the stability of black branes (black holes) in third order 
Lovelock gravity in that the stability can constrain the higher derivative coupling constants. 
The instability of D-dimensional asymptotically flat Einstein-Gauss-Bonnet and Lovelock 
black holes has been discussed by several authors [35-38]. Their results show that for 
gravitational perturbations of Schwarzschild black holes in D > 5 Gauss-Bonnet gravity, 
instability occurs only for D = 5 and D = 6 cases at large value of a' [36]. In [37], the 
authors showed that small black holes in Lovelock gravity are unstable. In this paper, 
we extend our previous computation to third order Lovelock gravity in I?-dimensional 
spacetime and show how stability constrains the Lovelock coupling constant. 

The plan of this paper is as follows. In section 2, we briefly review the thermodynamic 
properties of Reissner-Nordstrom-AdS black brane solution in third order Lovelock gravity. 
In section 3, we compute the viscosity to entropy density ratio via Kubo formula and its 
charge dependence. In section 4, the causahty problem is discussed. We study the stabihty 
issue of Reissner-Nordstrom-AdS black branes in third order Lovelock gravity in section 5. 
Conclusions and discussions are presented in the last section. 

2 Reissner-Nordstrom-AdS black brane in third order 
Lovelock gravity 

We start by introducing the following action in D dimensions which includes Lovelock 
terms and U{1) gauge field: 

/ = J d^xv^ ( - 2A + £i + + a'sjCs - AttGdF^^F'^^^ , (2.1) 

where 

Li = R, 

jC2 = R^.ysR''''^' - ^R^uR""" + R\ 

^24i?M--«i? + IQR>'-R,,R''^ - URR^'^R^, + i^^ (2.2) 

A is the cosmo logical constant, ^3 Gauss-Bonnet and third order Lovelock 

coefficients, respectively. The field strength is defined as Ff^^{x) = df^A^{x) — duA^lx). 

3 



The thermodynamics and geometric properties of black objects in Lovelock gravity were 
studied in several papers [39-42]. From the action (12. ip . we can write down the equation 
of motion [43], 

E ■ ■ ■ Kit = s-^^^t, (2.3) 

k=0 

where T'^ = Fj^Fy^g^"^ — ^g^^^Fp^FP" . Note that for third order Lovelock gravity, we must 
deal with D- dimensional spacetimes with D > 7. 
If we choose 

.,2 



a 



a 



an 



(D-3)(Z^-4)' 3(D-3)---(D-6)' 
the charged black hole solution in D dimensions for this action is described by [38] 



(2.4) 



ds^ = -H{r)N^dt'^ + H-\r)dr^ + —hijdx'dx^, 

t 

Q 



A, 



47r(D-3)r^-3' 



(2.5a) 
(2.5b) 



with 



H{r) = ko + -<l 
a 



2;2 



+ 



qH 



1/3- 



A 



(D-l)(D-2) 
2P 



where the parameter / corresponds to AdS radius. The constant A^^ will be fixed later. 
Note that the constant value of can be ±1 or and hijdx^dx^ represents the line element 
of a (D — 2)-dimensional hypersurface with constant curvature {D — 2){D — 3)kQ and volume 
Vd-2- The gravitational mass M and the charge Q are expressed as 

{D - 2)Vd-2 



M 

Q 



IGttG 



-m, 



D 



2_ 27r(D-2)(D_-3)^2 



G 



D 



Taken the limit a'g, a'^ with k^ = 0, the solution corresponds to one for Reissner- 
Nordstrom-AdS (RN-AdS). The hydrodynamic analysis in this background has been done 
in [44,45]. 

One may notice that here we use a black hole solution by choosing particular values 
of a2 and a'^ so that our computation can be simplified greatly. Eq.f l2.3p with the choice 
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(12. 4p yields one real and two complex solutions. We use the real solution in fl2.5ap . The 
general solution of third order Lovelock gravity in D dimensions for any arbitrary values 
of a'2 and ^3 was obtained in [38], but the line element of the metric turns out to be very 
complicated. Furthermore, the general solution may present naked singularities, which is 
not what we are interested [38]. In this paper, we only focus on the special case given in 

In the following, we mainly focus on dimensional case with fco = 0. Defining A = a/Z^, 
the function Hir) becomes 



j.D-1 ^D-l ^2D-A 

1 - 3A ( 1 - - a-±j-^ + a-±jj^^ 



1/3- 



:2.6) 



where a = jj-i ■ The event horizon is located at r = r+. The constant A^^ in the metric 



fl2.5ap can be fixed at the boundary whose geometry would reduce to flat Minkowski metric 
conformaly, i.e. ds^ oc — c^dt^ + dx"^. On the boundary r — > 00, we have 

so that A^^ is found to be 



N'^ = (2 7) 

l-(l-3A)i/3- ^^-'^ 

Note that the boundary speed of light is specified to be unity c = 1. From (12.61) . one 
can assume A < 1/3 in order to work with the branch of well-behaved solutions. That is 
because when A > 1/3 the sign of H{r) becomes minus in the asymptotic infinity, and we 
cannot recover the AdS geometry [*]. In section 4, we will carry out the causality analysis 
and find the causality constraints imposed on the value of A. 

We shall give thermodynamic quantities of this background. The temperature at the 
event horizon is defined as 



2TTy/g^ dr AnP 



T = ^ Y: = TZyi - 1) - P - 3)a] • (2. 



*In the Gauss-Bonnet case, the function H{r) has a different from: H{r) — 

1 1 ^ f — 4Agb(1 — 7t)(l ~ 7t)(1 ~ 72) j" I which imphes that the significant value of A lies 

in the range Ags < 1/4. Beyond this point, the Einstein-Maxwell-Gauss-Bonnet action does not admit 
a vacuum AdS solution, and then the AdS/CFT correspondence is undefined. In [27], it was found that 
causality requires exactly Xgb < 1/4 in the _D — s- 00 limit. This result matches precisely the assumption 
(i.e. Xgb < 1/4) used in [16,32]. 



5 



The black brane approaches extremal when a — * (i.e. T — 0). The entropy of black 
branes with ko = obeys the area law [38] and thus the entropy density has the form, 

3 Viscosity to entropy density ratio 

We explored the charge dependence of rj/s in the presence of Gauss-Bonnet terms for D- 
dimensional AdS black branes in [27]. In this section, we generalize the previous result on 
rj/s [26,27] to third order Lovelock gravity. It is convenient to introduce coordinate in the 
following computation 

r P _ f - 

z = uj = —uj, ks = —ks, f{z) = —H{r), 
r_|_ 

m = l^fl - (l - 3A(i - ^ + ^))'''] (3.1) 



2A 

We now study the tensor type perturbation h^^(t,X3, z) = (f)(t,X3,z) on the black brane 
background of the form 

/ D-2 
I Orhd Tr, r mTi rl'Tr, -I- 

/(^) P 



ds'^ = -f{z)N^dt'^ + —— + -- 20(t, X3, z)dxidx2 + ^ dx] 



i=l 



Using Fourier decomposition 



and expanding the action for tenor type gravitational perturbations (j){t, X3, z) to the second 
order, we obtain the effective action in the momentum space 

1 f dujdks 



^ = 16^ J (2^0^^^^ (M(^)0V + M,{z)<p') , (3.2) 

where the prime denotes the derivative with respect to z. 

An easy way to obtain the equation of motion of the tensor type perturbation is to 
substitute the fluctuated metric into Eq. (12.31) . One then find the hnearized equation of 
motion for (f){z) from the third order Lovelock field equation: 

M{z)(j)"{z) + M'{z)(l)'{z) + M2(f){z) = (3.3) 
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where 



/ {l - ^ [--V + z-\D - 5)/] + ^ [2f +{D- 7).-V'] /} 



2 



M2 = M{z)^^-klz^~'^ 



{l - p _ p _ 4) (/" + P - 5)P - 6).-V + 2p - 5).-V') 



+ 



2A 



2 



^-'(r + //") + 2(1^ - 7)^-V/' + \{D -7){D- 8)z-'f 



(3.4) 



[D-3)iD-A) 

We would like to emphasize that when D = 5 and the terms vanished, fl3.3p reduces 
to the main equation obtained in [16,26]. The shear viscosity involves physics in the 
lower frequency and lower momentum limit and one can neglect the M2{z) term in solving 
Eq. (13.31) . For the convenient calculation of the shear viscosity, we would like to introduce 
a new variable u = ^ and rewrite equation (13.31) in the new coordinate 

J{u)(f)"{u) + J\u)^'{u) + J2{u)<p{u) = 0, (3.5) 

where J{u) = and J2{u) = M2{\). 

In order to solve the equation of motion (13.51) in hydrodynamic regime, let us assume 
that the solution yields 

<P{u) = {l-uY F{u), (3.6) 

where F [u] is regular at the horizon, u = ±i^f can be fixed by substituting (13. 6p into 
the equation of motion, which we choose 

CO 

47rT' 

Since we only need to know the behavior at u ^ region, it is sufficient to expand the 
solution in terms of frequencies up to the linear order of u 

F{u) = Fo{u) + vF^{u) + C(z/^ kl). (3.7) 

The equation governing Fq{u) goes as 

[J{u)F^{u)\ = 0, (3.8) 

and can be solved as 

F^iu) = ^, (3.9) 



where Ci is an integration constant and must be zero as J{u) goes zero at the horizon so 
that Fo{u) is regular at the horizon. Therefore, Fo{u) is a constant, i.e. Fo{u) — C. Prom 
the equation at 0{v), 

[J{u)F[{u)\ - (i^'^h) = 0' (3-10) 



we find that the solution can be written as 



Regularity of Fi{v) at the horizon requires that 
C2 = - 



{{D -1)-{D- 3)a) (1 - -^{{D -1)-{D- 3)a)) C. (3.12) 



The value of C can fixed by the boundary condition C = lim^^o^l'i*) = 1- It is worth to 
noting that the above calculation is same as the Gauss-Bonnet cases given in [26,27]. 
Using the equation of motion, we write down the on-shell action 



1^ on— shell 



u=l 



The shear viscosity can be read off using the Kubo formula 

^. ImG(a;,0) r^~'N J{u)(j>{u)(l>'{u)\^=o 
77 = lim = — ^ lim 

D-2- 



(3.13) 



u=0 



^(7fcj)(^-ir^K^-i'-(^ (3.14) 
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The ratio of the shear viscosity to the entropy density turns out to be 

We obtain the same result as that of [27], which is also consistent with the prediction made 
in [19] when a = 0. In other words, the third order Lovelock coupling constant a'^ (or 
in our case) does not contribute to the shear viscosity. 

4 Causality 

The shear viscosity above is calculated in the hydrodynamical regime ( i.e. k^^ 0). In 
this and next sections, we will turn to the high frequency regime ( i.e. k'^ — >■ oo) and 
discuss the causality issue. 



Due to higher derivative terms in the gravity action, the equation (13.31) for the prop- 
agation of a transverse graviton differs from that of a minimally coupled massless scalar 
field propagating in the same background geometry. Writing the wave function as 

=e-*^*+^^^+*^3^^ (4.1) 

and taking large momenta limit ^ oo, one can find that the equation of motion (13. 3p 
reduces to 

fc'^r^^^ ^ 0, (4.2) 

where the effective metric is given by 



(4.3) 



Note that can be interpreted as the local speed of graviton: 



- ^'^ (4.4) 

where 



hi = {l- p _ 3)(^ _ ^) if" +{D-5){D- Q)z-'f + 2{D - 5)z-'f') 



+ 



2A 



2 



^-'{f" + ff") + 2(/^ - 7)z-'ff + 1{D -7)iD- 8)z-'f 



{D-3){D-A) [ ' ' ' 2 

= { 1 - ^ [z~'f + z-\D - 5)/] + ^ [2f + {D- l)z-^r\ /} . 

The above equations can exactly reduce to Gauss-Bonnet cases found in [16,17], if the 
terms vanished. For D = 10, we can expand near the boundary ^ = 0, 

2- 4[l-(l-3A)^/3]' 1 

^~ 4(l-3A)2/3 + 2(l-3A)V3_i +^^^^- ^^-^^ 

We can see that — 1 is always negative, which implies that the local speed of graviton is 
smaller than the local speed of light of the boundary CFT without any charge dependence 
. One can further check that for D = 7,8,9..., the first term in — 1 is also negative. 
Therefore, from (14. 5 p we can see that the causality imposes no constraints on possible 
values of A. Figure 1 demonstrates that the value of cl lies in the region 0.80 < < —0.8 
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as a function of u and A. Following the procedure of Ref. [17], one can find that the group 
velocity of the graviton is given by 



dk 



(4.6) 



So different from the Gauss-Bonnet cases [17,26], there is no causality violation in third 
order Lovelock theory with the particular choice of ctg and . The difference comes from 
the fact that terms change the causal structure of the boundary CFT. For more than 
third order Lovelock theory, the causal structure might be further modified by a[{i > 3). 




Figure 1: as a function of u and A when we choose D = 7 and a = 1.4. The hnes correspond 
to 0.8,0.6, —0.8, respectively, from left to right. 



One may notice that our discussions on the causal structure of third order Lovelock 
gravity are based on the metric fl2.5al) and the equation motion in high frequency limit 
(14. 2p . Hence, if we consider the general solution of third order Lovelock gravity with 
arbitrary values of a'2 and 0^3, we may find totally different causal structure. It remains to 
be carried out in the future. 



5 Instability 

In section 4, we have demonstrated that for the RN-AdS black brane in third order Lovelock 
theory, causality violation does not happen, which implies that the results obtained in 
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[17, 26, 27] might not be so universal as we expected. In this section, we extend our 
previous work on black brane stability to third order Lovelock gravity. 
Now, we rewrite the main equation in a Schrodinger form, 

__ + V' (,(..)) ^ = ^=^y (5.1) 

where ip {z{r^)) and the potential is defined by 



Vi{z) = 



dlnK(z)V . d f „, ^dhiKiz) 



dz ) dz \ dz 



(5.2) 



In the large momentum limit, the Schrodinger potential develops a negative gap near the 
horizon and the negative-valued potential in turn leads to instability of the black brane. 
In the large momenta limit A;'^ — > oo, the potential is mainly contributed by k'^c^g- For 
charged black branes, can be negative near the horizon and the potential is deep enough 
to have bound states living there. The negative-valued potential yields negative energy 
eigenvalue (i.e. u"^ < 0). The imaginary part of oj can then be positive. Substituting 
the eigenvalue of uj to the wave function for tensor type perturbations, one immediately 
find that perturbations grow as time goes on and the black branes thus are unstable. The 
negative-valued energy bound state corresponds to modes of tachyonic mass Minkowski 
slices and signals an instability of the black brane [46]. Let us expand in series of 
{l-u), 

c] ^N^ [(D -1)-{D- 3)a] {D^ [2X\a - if + 2{a + 1)A - l] 

-D [4A^(3a^ - 4a + 1) + 2A(a + 7) - 7] + [X'^(3a - if - 6A(a - 1) - 6] } 
{(D - 4) [{D - 3) - 2A (L» - 1 - (D - 3)a)]}-' {I - u) + 0{{1 - uf). (5.3) 

Since < a < and < w < 1, will be negative, if 

{D^ [2A2(a-l)2 + 2(a + l)A-l] 

-D [4A2(3a2 - 4a + 1) + 2A(a + 7) - 7] + [\\3a - if - 6A(a - 1) - 6] } 

{{D - 4) [{D -3) -2X{D -1- {D - 3)a)]y^ < 0. (5.4) 
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From the above formula, we find tlie critical value of A, 

X,= ^{-iD- 1){D -Q)-{D-3){D + 2)a 
+ |(L> - 1)2(3L>2 - 26D + 60) + {D- 3f{3D^ - lOD + 28) 

1 _2 

-2a(L)^-14D^ + 79D2-174D + 108)|'||L'-l-(D-3)a| . (5.5) 

Above the line of Ac, can be negative. Eg. (15. 5 p tells us that the stability of the black brane 
depends on the charge. The minimal value of Ac can be obtained in the limit a ("§z|); 

_l (D-3)(D-4) 

When D = 5, we recover the result found in [26]. Usually, for the application of AdS/CFT 
correspondence, we do not need to take infinite dimensionality limit. But the stability 
of higher dimensional black holes itself is an important topic in the study of black hole 
physics. 

The Einstein-Hilbert action is just the first term in the derivative expansion in a low 
energy effective theory. The Gauss-Bonnet and the third order Lovelock terms can be 
regarded as higher order corrections to the Einstein gravity. In this sense, the higher 
derivative gravity coupling constants should be small. In our discussions, we have found 
that the coupling constant \ {= {D — 3){D — 4)a//^) depends on the dimensionality D. 
But it seems that for fixed a and AdS radius / as D approaches infinity, A would be very 
large. That is not what we want. By doing stability analysis, we will find a way to restrict 
the value of A. As the value of D increases, one finds that Ac, min is bounded by 1/4, 

lim Ac = 7 (5.7) 

(i3,a)^(oo,g^) 4 

Thus we reproduce the result of [27]. Third order Lovelock gravity in our case does not 
add new constraints on the stability of the black brane. To show explicitly the behavior 
of gravitational perturbation in higher dimensions {D > 7), we solve the Schrodinger 
equation (15.11) with negative valued potential numerically and find unstable quasinormal 
modes (QNMs)(see tables [H [2] and [3] ). 

Table [U demonstrates that the unstable modes of the black branes are suppressed as D 
increases. Table 2 and 3 tell us the same story as we found in [27], that is to say, lower 
value of charge- (a) and A stabilize the perturbation, while the lower value of D strengthens 
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U 


A — u.oo 


A — u.ou 


A — U.ZD 


A — U.ZZ 


A — U.iO 


7 


oo. tc>6\j% 


o4.Z / (ill 


OQ CC/17^' 

ZO.0D4/Z 


91 Qni^7.)' 
Zi.oUO / z 


1 n nn i^fi.)' 

iU.UUOoZ 


8 


27.8612i 


31.3230i 


22.1646 


14.6227i 


5.1492i 


9 


23A12U 


22.7276i 


17.9269i 


11.2243i 


3.1228i 


10 


22.0265/ 


18.5886/ 


12.2811/ 







Table 1: Unstable QNMs for third order charged Lovelock black brane perturbation of tensor type 
for fixed charge (a = 1.20) and ks = 500. As D increases, the unstable modes are suppressed. 
And also, small A helps to smooth the perturbation. 



A 


a = 1.3 


a = 1.2 


a = 1.0 


a = 0.8 


a = 0.6 


0.33 


31.6732i 


27.8612i 


12.7312i 






0.30 


31.5500i 


27.3852i 


9.9506i 






0.27 


28.1982i 


23.7183i 


5.1296i 






0.24 


23.4156i 


18.6473i 









Table 2: Unstable QNMs for third order charged Lovelock black brane perturbation of tensor type 
for fixed dimensionality {D = 8) and ^3 = 500. This table indicates that instability is increased 
by a chemical potential. 



D 


a = 1.4 


a = 1.3 


a = 1.2 


a = 0.8 


7 


42.7952i 


39.3407i 


33.3140i 




8 


30.8934i 


29.5030i 


25.1211i 




9 




23.4155i 


20.6184i 




10 






17.442i 





Table 3: Unstable QNMs for third order charged Lovelock black brane perturbation of tensor type 
for fixed A (A = 0.28). This table shows combined effects of D and the chemical potential. Note 
that a = 1.4 exceeds the maximal value of charge permitted for 9- and 10-dimensional charged 
black brane and thus we leave the frequency blank there. 
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the instability. The reason for why higher D suppresses the gravitational fluctuation is 
because that no matter how big D is, A is bounded by 1/4, which means that for fixed 
AdS radius /, a' — > as the value of D goes up. The upper bound of A constrains the 
gravitational perturbation in the larger D limit. For QNMs of RN-AdS black holes in 
Einstein and Gauss-Bonnet gravity, one may refer to [36,47]. 

It would be very interesting to check for fixed value of charge, for which value of A 
the black brane becomes stable. In order to do this, one should first fix D in (15. 5p . then 
obtain a formula between A and a. Actually, (15.51) indicates that for A < Xc{D, a), the black 
brane becomes stable. For 5-dimensional black brane with charge in Gauss-Bonnet gravity, 
constraints from causality as well as stability separate the physics into four regions in (a. A) 
space: consistent region; only causality violation region; only unstable modes region; both 
causality violation and unstable modes region (see figure 4 in [26] for details). But for 
the particular case we are considering here, since causality violation does not occur, we 
have only two phases in the (a. A) space: stable and unstable modes regions marked by 
(15. 5p . One thing one need to be aware of is that instability of the black brane does not 
correspond to any fundamental pathology with the theory. This is quiet different from the 
causality violation which means that a theory is pathological. In the dual gravitational 
description, the unstable QNMs is identified with unstable uniform plasma with respect to 
certain non-uniform perturbation [34]. 

6 Conclusions and discussions 

In conclusion, we derive the main equation for tensor type perturbation in third order 
Lovelock theory and compute the shear viscosity. The result turns out to be in agreement 
with the prediction made in [19] when a = 0, that is to say, the third order Lovelock term 
does not add new ingredients into the shear viscosity of Gauss-Bonnet theory. 

We notice that an interesting point comes from the causality analysis. While in the 
Gauss-Bonnet theory, causality could be violated in the boundary CFT, we do not find 
causality violation in third order Lovelock theory. From (14. 4p . we can see that the local 
speed of graviton depends on both a'2 (~ A) and a'^ (~ A^). Although we are working only 
with a special choice of and ^3, Eq. (14.41) implies that causality receives corrections 
from the a'^ term. Thus, the causal structure in general third order Lovelock gravity must 
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be different from the Gauss-Bonnet gravity. We also expect that higher than fourth order 
Lovelock theory may impose more constraints on the causal structure of the boundary 
CFT. 

The instability of charged black brane with third order Lovelock theory shows the 
same properties as that of Gauss-Bonnet corrections. We find that higher D suppresses 
the unstable modes, but larger value of charge and A strengthen the perturbation. As D 
approaches infinity, the stability requires A to be bounded by 1/4. This is an important 
observation in that Eq. fl2.6p indicates that A could be as big as 1/3 without any causality 
violation happens in third order Lovelock gravity. But A ~ 1/3 violates the assumption 
Agb < 1/4 used in [16,32]. Fortunately, after imposing the stability constraint, we can 
recover the requirement A < 1/4 and thus third order Lovelock and Gauss-Bonnet gravity 
are somehow consistent. 
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